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SUMMARY 


An analytical method has been developed which is capable of es- 
timating the chamber and throat conditions in a nozzle with a low 
radius of curvature throat. The method has been programmed using 
standard FORTRAN IV language and includes chemical equilibrium cal- 
culation subprograms (modified NASA Lewis program CEC71) as an integral 
part. Program CEC71 is described in Reference 1. The method determines 
detailed and gross rocket motor characteristics in the presence of 
striated flows and gives detailed results for the motor chamber and 
throat plane with as many as 20 discrete zones. The method employs 
a simultaneous solution of the mass, momentum, and energy equations 
and allows propellant types, 0/F ratios, propellant distribution, 
nozzle geometry, and injection schemes to be varied so to predict 
spatial velocity, density, pressure, and other thermodynamic vari- 
able distributions in the chamber as well as the throat. Results for 
small radius of curvature (down to R c /Ry=.25) have shown good comparison 
to experimental results. Both gaseous and liquid injection may be 
considered with frozen or equilibrium flow calculations. For parametric 
studies it is also persible to elect frozen flow with calorically perfect 
gases in order to conserve computational time. In this case it is 
possible to study the effects of complete mixing compared to striated 
flow. 
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INTRODUCTION 


Objective 

The objective of this study has been to develop a computer program 
for determining the performance of rocket engines in the presence of 
striated flows and small radius of curvature throats and to perform 
parametric studies to determine the influence of engine parameters 
on mass flow rate, specific impulse, and chamber pressure. The pro- 
gram has been equipped with subprograms (a modified form of CEC71 
computer program of Reference 1) to determine the thermochemical 
composition for steady, inviscid flow in discrete zones. The program 
has the capability of analyzing the chamber flow field for specified 
injector schemes with gaseous or liquid propellant injection, and of 
analyzing the inviscid, transonic throat plane flow field. Predictions 
are made of the chamber operating conditions and the spatial distri- 
bution of velocity and thermodynamic variables for the reactants and 
products of combustion. Parameters which are considered for studies 
include the propellant types (especially gaseous 0 2 and H 2 ), 0/F ratios, 
propellant mass distributions, nozzle geometry, and injection schemes. 
Approximate two dimensional performance parameters at the throat plane 
are determined and compaired to one dimensional performance results. 

The philosophy used to develop this computer program is discussed in 
Appendix A. 



Considerations and Assumptions 


In small rocket motors heat transfer can be significant in terms of 
engine efficiency as well as thermal protection requirements. One method 
of reducing heat transfer is to reduce the heat transfer coefficient. 

This can sometimes be accomplished by laminarization of the boundary lay- 
er by using a highly convergent subsonic section of a nozzle. Experimental 

2 3 

work * has shown that there is a tendency with a rapidly converging sub- 
sonic section to stablize the boundary layer as a result of the rapid ac- 
celeration of the working fluid. This tendency toward laminarization of 
the boundary layer through the transonic throat region may be continued 
by employing a low throat radius of curvature ratio throat. In addition 
to the effect of reducing the heat transfer a decrease in nozzle surface 
area can be achieved. Further reductions in overall heat transfer to 
the nozzle walls may be effected by use of barrier cooling. By surrounding 
the hot core flow with fuel or ozidizer rich flow, temperature gradients 
to the motor walls are reduced thus the heat transfer is reduced. Stri- 
ating therefore permits maximum temperature in the core zone and yet 
permits near maximum overall specific impuse. 

In the treatment of striated flow with throat radius of curvature 
effects no closed form relationship exists between the chamber pressure 
and the mass flow. This is due to variations of total enthalpy and 
velocity in the chamber because of striation, as well as, radial pressure 
gradients at the physical throat due to radius of curvature effects. 

Thus choking with a given throat radius can only be assured and the flow 
parameters determined with a coupled solution of the throat and the chamber. 
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Even for unstriated nozzle flow the solution for the two dimensional 

flow field with small radius of curvature ratio throats is difficult and 
4 

time consuming . With inclusion of striation the problem becomes more 
complex and determination of performance parameters for a parametric 
analysis may require considerable machine computations. Radial velocities 
at the physical throat plane are thus assumed to be small and the throat 
is a revolved circular arc section. Figures la and lb present the flow 
regime of interest in which mass, momentum, and energy are conserved along 
streamlines. This analysis is for inviscid, compressible, axially sym- 
metric flow. Each annular flow zone is assumed to flow isentropically 
after leaving the chamber and satisfy conservation equations of, mass and 
momentum at the throat. Losses are limited to curvature and striation 
effects and to the degree of equilibrium flow attained (kinetic or energy 
conversion losses). The energy conversion losses are dependent on the 
degree of freezing considered for the solution and is restricted to lim- 
iting cases. In this analysis it was not necessary to average the mass 
flow nor to require the isentropic coefficient to be constant in a flow 
zone or stream tube at the throat. 

The chamber static pressure is the pressure at which chamber combustion 
occurs. This pressure is used for thermochemical calculations for each 
flow zone. In calculating an individual flow zone total temperature the 
difference between the kinetic energies of the burned gas and the injected 
reactants in the chamber zone is not considered since the difference is 
small compared to the overall energy of the material in the combustion zone. 
Thus, even for low contraction ratios the stagnation temperature in a flow 
zone determined with the chamber static pressure differs very little from 
the exact chamber total temperature for the zone. 
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THEORY 


1.0 Chamber Solutions 


The chamber solution depends on the method of mass generation and 
on the properties of the injected and combustion materials. To obtain 
choked conditions at the throat one may have to alter the chamber combus- 
tion static pressure or the total mass flow and recalcualte throat plane 
parameters. By this technique the coupling is maintained between the 
chamber solution and the throat plane solution. By multiple iteration 
the complete solution can be made to simultaneously satisfy chamber 
and throat requirements. 

For a discrete number of stream tube zones as shown in Figure la, 
the mass, momentum, and energy equations for one dimensional, inviscid, 
compressible flow are: 


Pt Wt + pi vv t A, nil + nii m. 

TU i A fu. OX i A OX. 1 OX i 1 fu i A OX.j 


( 1 . 1 ) 



m T w T + m T Wt + A P T 
1 fu i 1 fu i 1 ox i 1 ox i 9 i 


P A 

g 


m . w 


where 



- h + w 2/2 
9 t 9i 




O 

dT + h 9i 


( 1 . 2 ) 

(1.3) 
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In addition, it is necessary that 



m T = (0/F) • m T 

I fu i 


(1.4) 

(1.5) 


Further, if the injected propellants are gaseous and perfect gases are 
assumed 


and 


fu i 


* p iA r i 


fu. 


fU i 



( 1 . 6 ) 

(1.7) 


For a constant chamber static pressure , the propellant distributtion 
in each zone is specified as the fraction of the total mass flow. This 
case is of interest when it is necessary to compare engine performance at 
a fixed pressure. The total mass flow is determined by iteration and is not 
known apriori. An initial estimation of the total mass flow or an adjust- 
ment by iteration of the total mass flow permits intermediate calculations 
and allows determination of parameters in the chamber. Thus, iteration 
between the chamber solution and the throat plane solution drives the 
intermediate results to the solution values of parameters. An analysis with 
a constant chamber static pressure is a practical mode to consider for 
parametric studies since the chamber static pressure, propellant reactants 
inputs, and 0/F in each zone, once specified, do not change results of 
chemical equilibrium calculations for the mass flow adjustments required 
to determine the solution. The solution being that for the chamber and 
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the throat plane with a specified throat radius of curvature. Equations 
(1.1) through (1.4) provide a system of equations to which the unknowns 
T , W , A and P. may be determined. 

9-j 9-j 9j 1 

For specified mass flow in each zone , there exists a chamber static 
pressure which simultaneously satisfies Equations (1.1) through (1.4). 

By application of the governing equations at the throat plane, the throat 
radius determined is the prescribed throat radius. For this scheme of mass 
generation, chemical equilibrium calculations would be required at each 
iteration step. Initial estimates of thermodynamic variables are altered 
once to provide data for intermediate calculations and a final solution. 

For specified injector reservior pressures , injection equations for 
mass generation due to the reservior pressure for the fuel mass flow injected 
into each zone is determined by 


m. 


= C 


fu i 


D I. 1 


[2 P j (P c 


fu i 


fu. 


fu i 


-P r )] 


1/2 


and thus 

= m T (1 + 0/F.) 
fu- 1 

l 

and 


m. 


= m. - m n 


ox. 


1 fu . 


( 1 . 8 ) 


(1.9) 


( 1 . 10 ) 


This scheme is limited in the program by the requirement that injector fuel 
areas are required to be specified. Intermediate results are determined 
by a double iteration since the injector pressure, Pj, is also determined 
within the chamber solution iteration loop. Equation (1.8) is not strictly 
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applicable to gaseous injection but is operational for liquid or gas in- 
jection in the computer program . 

For cases where frozen, calorically perfect flow is considered, a solu- 
tion for fully mixed flow may be included. The results give a limiting 
case for comparison to striated flow. Appropriate thermodynamic and flow 
variables may be calculated as follows: 


• • 


C p = * C P V m total 


’mix 3 1 

m. 


MW. 


1 / 1 ^ • 


1 


"•total MW 9i 


R n = <$U/MW 


"mix 


c ft 


'mix 


Y g . ■ C P - R g . > 

ralx W 9 m1x 9|mx 


r <, ' 1 Cr „ T o„ *1 

9 mix 9 i 9 i 


C P "'total 


3 mix 


Conservation of energy for this case requires 


C P T o = h g . + w g . / 2 

9 mix 9 mix 9mlx 9mix 


but 


h g . “ C P T g . “ C P P g . 

mlx 9 mix 9mlx 9 mix-^2x 

p 9 . R g, 


mix 


mix 


= C D P n W n A n 

‘ g mix J jn i x g mi , x 9 
" total Rg mix 


( 1 . 11 ) 

( 1 . 12 ) 

(1.13) 

(1.14) 

(1.15) 


(1.16) 


(1.17) 
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so 


C D T = C D P n v/ A + w n /2 
P o P q . q . q q . 

g„. g . q„. ^mix y mix y y mix 

3 mix *mix 3 mix 


rti. . al R 
total g 


mix 


0.18) 


Conservation of momentum requires 


P A +m. + n w = P A + I m. w 
g . g total g . g g 1 q. 
a mix s y mix 3 3 y i 


(1.19) 


or 


"mix 


I m i w 

P + 9 i 

g —5 — 


m* . , w 
total g . 

y mix 


( 1 . 20 ) 


Elimination of P from Equation (1.18) with Equation (1.20) gives a 
9 mix 

quadratic in w . Equation (1.20) gives P and p may be 
9 mix 9 mix 9 mix 


determined by 


A1 so , 



m total^g w g 


mix 



( 1 . 21 ) 


( 1 . 22 ) 


If only one zone of gas flow for a low temperature or nonreacting gas 
is the subject of investigation, a simplified calculation option in the 
program may be used. Injector details are not determined and thermochemical 
data are input as opposed to determined by chemical equilibrium calculations. 
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The energy equation, mass equation, and equation of state are simultaneously 
solved. By elimination of the injector pressure, Pj, from the list of 
unknowns, the momentum equation is not required. This option in the 
program serves to predict radius of curvature effects at the throat plane 
for low temperature single zone gas flow. 

2.0 Throat Solution 

The properties of the combustor flow field are determined with the 
combustion chamber calculations prior to the throat solution. Initially, 
estimates are made followed by iterative corrections until the converged 
solution results for desired fixed parameters. The momentum and mass equa- 
tions describe the flow at the throat. These equations together with the 
requirement that energy be conserved in stream tubes defines the flow. 

A schematic representation of the geometry of a typical nozzle throat is 
given in Figure lb. 

The radial momentum equation is: 

3P/3r = - p[u 3u/3r + w 3u/3z - v^/r] (2.1) 

Considering the throat plane, we expect that the radial velocity will 
be small everywhere since it must be zero at the throat wall and at the 
nozzle centerline. Thus, for axially symmetric flow 

dP/dr = -pw du/dz (2.2) 

Since dr/dz = u/w, then on the throat plane between the circular arc 
throat wall and the centerline 
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dP/dr = -pw 2 d 2 r/dz 2 = -pw 2 /& 


(2.3) 


In the above equation Prefers to the radius of curvature of the stream- 
lines at the throat. The energy equation applied at the throat where 
u = 0 requires 

w 2 = 2(H o -h) (2.4) 

For frozen, calorically perfect flow, the isentropic coefficient, y, 
is constant in each stream tube and it is convenient to express the enthalpy 
in stream tubes as: 

l/Y 

h=^rr[P] y [c] (2.5) 

Where the constant C for each stream tube is defined by the isentropic 
process equation as: 

C = CU) = P/Jpf (2.6) 

The equation for mass conservation is: 

d*/dr = d^/dr (l/m tota ,)- (2.7) 

total 

Applying Equation (2.4) to Equations (2.3) and (2.7) it can be seen 
that a set of two nonlinear total differential equations result. The solu- 
tion of these equations can be obtained at the throat providing C, Ho, y, (V 
can be specified as a parametric function of \p. Of course the variation of 
C, H q , and y with ip is known for the current iteration from the chamber 
solution (at least for frozen, calorically perfect flow: for equili- 

brium or frozen flow these parameters are determined as a function of 
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pressure for a given chamber solution). It remains only to specify a re- 
lationship for R/ to close the formulation at the throat plane. It is clear 
from experimental and analytical studies that the main functional form for 
should be R/= R/(R C > R j,r) . This formulation for R/ neglects effects due to 
inlet geometry between the chamber and the throat region; however, experi- 
mental evidience indicates this to be a minor variation except for contrac- 
tion section wall slopes near perpendicular to the nozzle axis and/or 
5 6 7 

R /R = 0. ’ ’ • The appropriate boundary conditions for R/ regardless of 

0 i 

striation effects are: 

R/ = R c at r = Rj 

R/ = ® at r = 0 


A plausible form for which satisfies these conditions is 

- M 1 + V>~ r 


R/ = R 


r 


( 2 . 8 ) 


Figure 2 presents a number of calculations for with respect to R c /Rj 
for various values of n. Based on these results a value of n = 3/2 is 
indicated. The representation chosen for R/can be further evaluated by 
comparing results for the variation of throat plane Mach number with radius. 
This will be presented later in the report. 

Equations (2.3) and (2.7) thus comprise the necessary equations, ex- 
pressed in a more general form than previously considered., since H q = H q ( 4 >), 
h = h(P, i|>), p = p(P, ♦), and w = w(H Q ,h) = w( p , if;). In general, H q U) 
is determined from chamber conditions and is calculated before equations 
at the throat are solved. If required, C(<ji) is also determined from chamber 
conditions. Throat plane variables are determined so that the throat radius 
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which just chokes the flow must be the prescribed radius. An iteration 
between the chamber solution and the throat plane solution allows starting 
estimations to be corrected to assure the final iteration predicts choked 
flow for the desired geometry and fixed parameters. For simultaneous solu- 
tion of the throat plane equations, Hamming's prediction-corrector tech- 
g 

nique was employed. This method proved to be efficient and stable and is 
initiated with the Range-Kutta method so to be self starting. 

3.0 Performance Solution 

Two types of performance calculations are made. The first type includes 
the effects of throat radius of curvature and striation. The total mass 
flow is given by 

■"total =In, i < 3 -" 

The throat vacuum specific impulse is determined by integrating the 

velocity and pressure profiles across the throat plane so that 

r=Ry r=Ry 

I . f wd* ♦ /PdA/m tota , (3.2) 

vac r=0 r=0 

The characteristic velocity is 

A T P o /m total} 

9-: 

c* = , (3.3) 

m total 

Within the program, intergrations are actually represented as finite sum- 
mations based on the fraction of the mass or area over which the required 
thermodynamic variable is representative. In Equation (3.2) the velocity 
at a point on the throat plane is assigned a mass fraction determined by 
adjacent values of ip and pressure is assigned the area over which it is 
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representative. Small independent variable step sizes input to the program 
permit the integration approximation technique used to be accurate. 

For comparison, the second type of performance calculations are based 
on the ideal or optimum distribution of Mach number uniformly equal to 
one on the throat plane. Although this distribution does not in general 

g 

occur for one dimensional striated flow , it serves as a limit to which 
striated flow with radius of curvature effects may be compaired. If frozen, 
calorically perfect flow is considered, the one dimensional mass flow is: 


m, 


total 


= I 


1-D 


★ 

A i P o 


9 i 


(r g. /T o 
1 9 


, 1/2 


1 


(3.4) 


where 



(3.5) 


If frozen or equilibrium flow is considered, the one dimensional mass flow 
is determined by 

m total = 1 A i* P g /c *l-D i 


where c*-]_ D is determined by the chemical equilibrium calculations and 
^ * 

is based on P . A i in either case is determined by the throat solution 
and is the area of a particular flow zone on the throat plane. In a similar 
manner the one dimensional vacuum specific impulse at the throat is 


i T o -i 

SPyac l-D '''total'- 9 ’ 9 ' < T o 

Si 

I P 0 (P./P- ) A* 

+ ° 9i 1 ° 9j 1-D 1 


m 


total 


1-D 


(3.7) 


13 



where 


1 + v. 


< T o_y T i>i- D - 


(3.8) 


(P i /P o >1-0 * P 0„ KlV 


(3.9) 


for frozen, calorically perfect flow. For frozen or equilbirum flow 


: SP =ZI sp, „ V m total 

vac i_[) vac 


(3.10) 


where I 1 . 


is determined by the chemical equilibrium claculations. 


The characteristic velocity for one dimensional calculations is 

C, -° •total, „ 


(3.11) 


4.0 Chamber and Throat Iteration Techniques 


A chamber solution iteration is necessary to assure that Pj, , 

T , and w satisfy conservation equations. The injector pressure, Pj, 

9i 9^ 

converges after an initial estimate is made by the program. The iteration 


equation employed is as follows: 


T A 2 
* 9, . 


P I • p o + (p r - V ‘-TC 1 ) 


(4.1) 


A double iteration was employed to determine T , A , and w . This tech- 

9 n - 9j 9j 

nique employes a simultaneous solution of the energy, mass, and momen- 
tum equations with the assumed value of P as initially estimated or, 
after the first chamber solution iteration loop is completed, as determined by 
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Equation (4.1). Convergence is achieved after the above procedure is 
complete and the following is determined to be true 



The injector pressure, Pj, is the pressure at each injector surface in 
the chamber and is the same for each zone. Within the program Pj is a 
dimensioned variable and is identical with the variable PIAVG. If conver- 
gence fails for some reason, variables P , PIAVG, 2A , and A c are printed 
in an error message. Convergence difficults may be encountered for a single 
zone of flow of low temperature or non reacting gas. (This is because w g 
becomes large for this case). Special subroutine SPPG has been provided 
in the program to avoid convergence difficulties should they occur for the 
above mentioned single zone flow case. 

To initiate the throat plane solution an estimate of the centerline 
pressure is required. The estimation is based on a small perturbation 
method^ 0 



With the initial estimate, throat equations are integrated by Hamming's 
g 

method to determine the throat wall radius corresponding to <p = 1.0. 

The pressure, P CL , given by Equation (4.3) is increased and a new wall 
radius is determined by integration. Following this, dr^dP^ is used to alter 
P^ L until the minimum radius is determined. The iteration process by which 
dr w /dPcL = 0 is sought requires that dP^ be sufficiently small and that 
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r w = R . Only seven iterations are permitted by the program once r 

is bracketed. A convergence message is provided by the program if seven 
iterations are not sufficient; however, calculations are not halted but pro- 
cede to alter or m^^ such that r w = Ry. The minimum wall radius 
as determined by the above procedure is used to alter the chamber static 
pressure by 


P g - P g< r w . 'V 
y y min 


(4.4) 


for fixed mass flow or fixed reservior pressure solutions. For fixed chamber 
static pressure, the total mass flow is altered to obtain r w = Ry by 


min 


m 


total = m total ^VV, * 


(4.5) 


min 


If attempts to determine P^ or m ^ :0 | :a i f° r the solution exceed seven, another 
convergence message is provided and parameters as described in the program 
output are, nevertheless, determined. Results may be erroneous but 
sufficient for parametric studies. The most flagrant error occurs if the 
P^ estimation is very inaccurate. This error may occur for R c /Ry £ -25 
so are seldom encountered. This usually causes the pressure, P, to become 
zero and internally generated messages occur and machine generated mes- 
sages may occur. Calculations are not necessarily halted and corrections 
provided in the program may allow the solution to be determined. 
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5.0 Results and Conclusions 


In this section equations which were presented under Throat Solution and 
Performance Solution will be more thoroughly qualified. Results of the present 
analysis are discussed and where possible compared with experimental data and 
with other methods of analysis. The analytical form used for the flow radius 
of curvature at the throat plane is presented and verified by comparison to 
experimental data. Results of the many options available in the program will 
not all be presented; however, several typical results will be presented and 
discussed. Results of parametric studies are presented as well as performance 
calculation results. 


5.1 Effect of Radial Velocity at the Throat Plane 


Consider an axisymmetric converging-diverging nozzle shown in Figure 2. 
An element of mass flow, dm, passes through surface dl inclined to the axis 
at angle <j>. The element of mass has velocity V inclined at angle a to the 
nozzle axis and is da wide. The elemental mass flow is 


Since 


dm = 2irrpVda 


da = sin(<)>-a) dl 

and dl = dr/si n<j> 

then, 

dm = 2 *p v sin U-g) rdr 
sin<}> 

On the throat plane. $ = tt/2 so that 


m 


total 


r=R T 

-f 


2irp Vcoso rdr 


r=0 


(5.1) 

(5.2) 

(5.3) 

(5.4) 


(5.5) 
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Thus, the discharge coefficient is 



(5.6) 


The effect of radial velocity on mass flow at the throat plane is exhib- 
ited in the density as shown by Equation (5.5). For single flow zone of a 
frozen, calorically perfect gas the density may be written as 


p 



vl 


2 2 
w +u 


<Wo, 



(5.7) 


The radial velocity component, u , may be assumed to be negligible at the 

2 

throat plane compared to w since u = 0 at r = 0 and at r = Ry. Using two di- 
mensional methods of Reference 4 this assumption may be verified for nonstriated 
flow for values of R^/Ry down to .35. It is therefore concluded that this as- 
sumption will be useful for the present analysis. Recall, however, au/3z f 0 
on the throat plane. 


5.2 Streamline Radius of Curvature at the Throat Plane 


In the discussion of the throat plane solution the radius of curvature 
of the streamlines was given by Equation (2.8) as 

. rV 1+ V r > - r l n 

R C J (5.8) 

This form satisfies the boundary conditions at r= 0 and at r = Ry for 
any n greater than zero. Figure 3 presents calculation results for n = 1.0, 1.5, 

and 2.0. Based on these results a value of n = 1.5 is indicated as the appropriate 
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exponent. In addition, Figure 4 shows good agreement between the present 
theory throat plane Mach number distribution and experimental data of Reference 5 
for R /R t = ,625. Other theories are also shown for comparison. Figure 5 
presents the distribution of throat plane static pressure to chamber pressure 
ratio of the present method. Results compare favorably to the data taken by 

5 

Cuffel , Back, and Massier. 

5.3 Results of Analytical Method 

The effects of striating mixtures of gaseous hydrogen and oxygen for 
nozzles with small throat radius of curvature to throat radius ratio having 
various amounts of individual flow zone mass flow have been investigated. 

Figure 6 illustrates the dependence of the mass discharge coefficient with 
respect to the cooling mass flow ratio for R c /Rj = .5 and P g = 200 psia. 

Results are for equilibrium; frozen; frozen, calorically perfect; and frozen, 
calorically perfect fully mixed flows. Figure 7 shows similar results for 
R c /Ry = 2.0. Figures 6 and 7 indicate that if striation and equilibrium are 
maintained, an improved discharge coefficient may be obtained. Improvement in 
discharge coefficient results for equilibrium flow compared to fully mixed flow 
expecially for a large mass of fuel rich cooling. Both Figures 6 and 7 show 
that little change in discharge coefficient with fraction of cooling mass flow 
occurs when chamber flow zones burned gases are fully mixed in the chamber. 

Figure 8 presents discharge coefficient, Cp, versus radius of curvature ratio, 
R c /R t , for various amounts of cooling with unmixed and fully mixed flows. Input 
thermochemical data used to calculate results shown in Figure 8 for frozen, 
calorically perfect flow at the throat plane were determined at 50 psia chamber 
pressure with constant chamber total enthalpy in flow zones. The core and cooling 
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zones had 0/F ratios of 6.0 and 0.5, respectively. As in Figures 6 and 7 striated 
flow has higher discharge coefficient than uncooled or fully mixed flow for these 
propellants. 

Figure 9 shows the spatial distribution of Mach number and static pressure 
on the throat plane for R c /Ry = 0.5, 1.0, 2.0, and 10.0 with four flow zones. 

The Mach number discontinuities between flow zones appears for all R c /Ry resulting 
in a discharge coefficient less than 1.0. A discharge coefficient of unity may 
be obtained only for a Mach number at the throat plane uniformly equal to 1.0. 

This point will be expanded later in this section. It also appears that striation 
may improve the discharge coefficient. This is verified by the hydrogen and oxygen 
propellant combinations used to determine Figures 6, 7, and 8. 

5.4 Effect of Throat Plane Mach Number Distribution on Discharge Coefficient 

Based on Equation (2.3) the axial velocity may be found by 

W 2 =-Bf dP/dr (5.9) 

Since u = 0, the density is p = p (P ,ip) and thus mass flow rate given by Equation 
(5.5) and discharge coefficient given by Equation (5.6) depend on pressure and 
the radial pressure gradient on the throat plane, on thermochemical parameters 
of chamber flow zone gases (on t|>), and on the flow radius of curvature, R/. For 
frozen, calorically perfect flow, the mass flow with radius of curvature and 
striation effects by the present analysis is 

; . & 7 t V7 rdr (5.io) 

tota js; J x/r 

v r=0 
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In the present analysis the one dimensional mass flow for frozen, calorically 
perfect flow was determined by. 


total 


1-D 


*R T 2 /V& 



WCTT-lT 

a i 

(5.11) 


Thus, by taking differential elements of radius with associated thermodynamic 
variables the discharge coefficient may be written as 






j 


M. 

J 



(5.12) 


For unstraited flow considerable simplification results and the discharge 


coefficient for frozen, calorically perfect isentropic flow is 


r - 2 

c °' 7 


? T 


V 1 


Y g +1 


+1 v STyT^T) 


(Vp, 


r=0 


r dr 
(5.13) 


In Equation (5..13) the coefficient is one if the Mach number is uniformly equal 
to one at the throat plane. For all other Mach number distributions consistent 
with conservation of mass, momentum, and energy the discharge coefficient is 
less than one. 
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That the discharge coefficient of striated flows in the presence of small 
radius of curvature throats tends to be greater than the discharge coefficient of 
unstriated flow is not necessarily true for all fuel-oxidizer combinations. 

Figure 10 presents R c /Rj versus discharge coefficient examples for comparison. 

In Figure 10 the hydrogen-oxygen propellant combination with a core 0/F = 7.0 
and a sheath 0/F = 0.0 shows an improved discharge coefficient for 20% striation 
over unstriated flow. However, for the hydrazine-nitrogen tetroxide combination, 
which also assumes frozen, calorically perfect flow, the unstriated flow dis- 
charge coefficient is greater than the striated flow discharge coefficient. As 
in Figure 8, Figure 10 employes thermochemical data determined for a chamber 
static pressure of 200 psia and constant total enthalpy. For Figures 8 and 10 
the isentropic exponent in the chamber core flow zone is not equal to the ratio 
of specific heats. The significant difference between the flows used to cal- 
culate results for Figure 10 is that for striated flow using and ^0^ 
propellant combinations less throat area is required to accommodate 
sheath zone mass flux than is required for the H 2 and 0 2 propellant combinations 
using H 2 as a coolant given the same fraction of coolant mass flow to total mass 
fl ow. 

5.5 Parametric Study Results 

Extensive parametric studies have consisted of determining the effect on 
discharge coefficient due to striation and radius of curvature ratio with frozen, 
calorically perfect flow. By the direct imput of thermochemical data and the 
use of a frozen, calorically perfect flow analysis, particular sets of thermo- 
dynamic variables were fixed and used to generate parametric curves with ther- 
modynamic parameter combinations. Ratios of specific heats, stagnation temperatures. 
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and molecular weights as well as cooling mass flow to total mass flow appear 
to be the main parameters required to develop a parametric relationship between 
striated and unstriated flow discharge coefficients for given throat radius 
of curvature ratio. Studies have yielded no general parameter combination which 
correlates striated to unstriated flow discharge coefficient without knowledge 
of the results of the throat plane solution. However, trends have been observed, 
namely, there is a significant effect of combustion gas ratios of specific heats 
in striated flow. Stagnation temperatures and molecular weights of the combus- 
tion gasses have only a secondary effect on discharge coefficient. Even for 
unstriated flow there is a difference between discharge coefficients with dif- 
ferent ratios of specific heats. Figure 11 shows that discharge coefficient 
increases with decreasing ratio of specific heats. In Figure 11 the molecular 
weight and total temperature of the burned gas remained constant for all curves 
of R c /Ry versus C D . 

5.6 Performance Results 

The characteristic velocity, c*, versus cooling to total mass flow ratio 
is shown in Figure 12. Curves are for fully mixed and unmixed flow for radius 
of curvature ratios, R c /Ry, equal to 0.5 and 2.0. Under the same flow mixing 
conditions, i.e. unmixed or fully mixed, and for the same fraction of cooling 
to total mass flow the characteristic velocity is greater for R c /Ry =0.5 
than for R c /Ry = 2.0. This is due to the decrease in total mass flow at speci- 
fied chamber pressure for R c /Ry = 0.5 compared to the total mass flow for 
R c /Ry =2.0. Recall that 

2 

= V (VWal ) < 5 - ,4 > 

m total 9 i ' ' 
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In Figure 12, fully mixed flow at a given R c /Ry has a characteristic velocity 
greater than unmixed flow for a fixed chamber pressure for the same reason even 
though total pressure decreases for fully mixed flow from unmixed flow. 

Figure 13 presents the characteristic velocity, c*, versus cooling to total mass 
flow ratio. Curves are for two zone striated and two zone, fully mixed flow 
with R c /Ry =2.0. Also shown is the resulting chamber static pressure, P g , 
which, for this case, was the same for the striated and mixed flow. Note that 
the static pressure for this engine drops off as the cooling flow is increased. 
Results for Figure 13 were obtained by calculations with constant reservior pres- 
sures and may not be compared on the same basis with results shown in Figure 12. 
Figure 12 and Figure 13 do agree, however, in indicating an improved character- 
istic velocity of mixed flow over striated flow. Results for curves of Figure 13 
show a decreasing total mass flow with increasing cooling to total mass flow ra- 
tio which agrees with a decreasing chamber static pressure. At the same mass 
ratio, the total mass flow for striated flow is greater than the total mass flow 
for fully mixed flow. 

5.7 Conclusions 

The theory presented here permits an extension to an analysis of one dimen- 
sional compressible striated nozzle flow. The throat plane solution gives an 
approximation to two dimensional axi symmetric flow and has successfully been 
used for analyzing the throat plane in converging-diverging type nozzles with 
striated flow for throat radius of curvature to throat radius ratios down to 
.25. The method has been programmed using standard FORTRAN IV and is fully 
documented in this report. The program is capable of analyzing the nozzle throat 
plane and rocket motor chamber with frozen, calorically perfect; equilibrium; 
or frozen striated flows as well as fully mixed frozen, calorically perfect 
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flow. By coupling the chamber and throat plane solutions, chamber and throat 
plane spatial dist. ributions of velocity and other thermodynamic parameters are 
uniquely determined for specified geometry. Several modes of injector description 
and mass generation are possible with the program. 

Results of the program agree favorably with experimental data and results 
of other computational methods, especially in predicting the mass discharge 
coefficient. No experimental results are presently available which permit distinct 
annular zone or striated flow to be compared to results obtained by this method. 

The dependence of discharge coefficient is primarily on the radial pressure 
gradient, the adiabatic or isentropic coefficient, and the amount of mass in each 
flow zone. Radial velocity is insignificant especially if streamlines on the 
throat plane have radii of curvature centered on the throat plane as proposed 
by the present analysis. 

It can be concluded that results obtained by this method may offer an 
advantage over other methods by being in agreement with experimental data and 
by being in agreement with but requiring less computational time than finite 
difference methods used to obtain the same data. Further, multi zone striated 
flow solutions may be obtained with the program in only seconds of computer time 
which permits rapid calculation of parametric performance calculations for 
candidate engine designs. 
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DESCRIPTION OF PROGRAM INPUT 


Program input information or data consists of six types. Three of the 
types are required for all problems and three are optional. The optional 
inputs allow two kinds of problems. The required types consist of the 
following: 

(a) Comment Cards (two cards) 

(b) Interger control card (program variable name LETOUT) 

(c) Namelist data which includes input parameters and option integers. 

The optional input data include (1) thermochemical data for frozen, calor- 
ically perfect flow problems input under nameslist EQU or (2 and 3) modified 
input data to the chemical equilibrium calculations with formats identical 
to the formats required for the corresponding in NASA program CEC71 J 
Table I and II are schematics of input cards containing namelist and other 
input information or variables required for the two kinds of problems which 
may be considered. Input format information may be seen in Appendix D and 
examples of input cards for several typical problems are given in Appendix F. 

Comment Cards 

Two comment cards are required for each problem case. Any information 
may appear on these cards anywhere in the eighty column fields. 

Control Card 

One integer control card is required for each problem case. The 
integer, read by an 15 format provides control for the following: 
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(a) An end to all computations 

(b) Printing of variables in namelist NUM 

(c) Supression of printing of namelist NUM variables 

(d) Stacking cases behind an initial case where thermochemical data 
determined by chemical equilibrium calculations does not change 
from the initial case. 

If the control integer is zero computations are ended and the program 
ends. For example, if the last three cards on a set of stacked data are 
blank, the program ends after the next to last data set (the last case) 
has been used for calculations. Three blank cards should be used to provide 
a normal exit since each case requires two comment cards. No further 
printed output results if the integer control card is blank or the integer 
control variable, LETOUT, is zero. If the integer control is greater than 
one but not equal to 71, an initial case is considered but variables in 
namelist NUM are not printed. If thermochemical calculations occured for 
a preceding initial case and these data will be used again in the present 
case, an integer control equal to 71 is input. If this is so (LETOUT 
equal 71) and thermochemical calculations have been made in a previous 
initial case, input data in THERMO and REACTANTS sets are not to be input 
again. If thermochemical data are to be input under namelist EQU and not 
calculated, LETOUT may be equal to 71 with the only effect being that 
parameters in namelist NUM will not be output under the namelist write. 

Namelist NUM 

The third type of input data is required for all cases. This data is 
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read under a namelist read input and includes the following: 

(a) Throat plane integration data (required by Subprogram HPCG) 

(b) Geometry data 

(c) Chamber operation mode data or mass generation data (three modes 
are possible) 

(d) Initial estimation data 

(e) Option control data 

(f) Injector data 

A more complete description of the above is given in Appendix D. 

Throat plane integration data are required by subroutine HPCG. These 
include the arrays or vectors PRMT and ERWT. PRMT are required for the 
first or an initial case where the throat radius is changed from a pre- 
ceding initial case. ERWT are optional and if not included in namelist 
NUM, ERWT are determined internally. PRMT may only be changed by input 
with namelist NUM where the dependent variable error weights, ERWT, unless 
specifically designated, are optimized for throat solution iterations. 
Generally, ERWT are not input. 

Geometry data include the throat radius (RT), the throat radius of 
curvature ratio (RCRT), and the chamber to throat area contraction ratio 
(EPSC). Unless changed, geometry data for stacked cases remain the same as 
initial case data. 

Chamber operation mode data consist of only one of the three following 

sets. 

(a) PG, PCTMO (PCTMzero) 

(b) MOT I, AGAMGU, AMWGU, ATQGU (ATzeroGU) 

(c) PRF, ESPG, CDINJ, PRO 
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The chamber static pressure and the fraction of total mass flow in 
each flow zone starting with the inner most zone in all cases comprise the 
first chamber mode data set listed above. The sum of all PCTMO must equal 
1 . 0 . 

When specified mass flow is considered, not only is the mass flow 
in each flow zone required but also initial estimations of the 
isentropic coefficient, the molecular weight, and the total temperature 
in each flow zone are required. Initial estimations are replaced by equil- 
ibrium calculations or data input under namelist EQU during the first 
chamber solution iteration. 

When specified reservior pressures are considered, the fuel and oxi- 
dizer reservior pressures are required. The 0/F ratio in each zone permits 
the oxidizer injector areas to therefore be determined. An initial esti- 
mation of chamber static pressure is required to start the chamber solu- 
tion and the throat solution iterations. The initial estimate of the cham- 
ber static pressure is determined by a fraction of the fuel reservior 
pressure in zone one. This estimate, ESPG, is optional in namelist NUM 
input but ESPG = .7 is set before data is read under namelist NUM. ESPG 
must be less than one and should be based on experience with rocket motor 
calculations to avoid extensive iteration calculations in the chamber and 
throat solutions. Suitable values of ESPG have been determined to be .9 
for gaseous injection and .5 for liquid injection to speed convergence some. 

. Option control data determine the type of problem solved and must be 
considered carefully to assure proper program operation. Integer variables 
as opposed to logical variables have been used to give more flexibility. 

Six integer variables are included in this category and are as follows: 
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MODECH 


(a) 

(b) NCEC71 

(c) NAORW 

(d) MGORLI 

(e) MTOTMX 

(f) NSPPG 

MODECH must be input for each initial case and need not be input again 
for stacked cases where it remains the same. This control allows specified 
mass flow for MODECH = 1, specified reservior pressures for MODECH = 2, 
and fixed chamber static pressure for MODECH = 3. NCEC71 has the value 0 
before read by namelist NUM. If NCEC71 = 0, frozen, calorically perfect 
thermochemical data is read by namelist EQU. NCEC71 = 1 permits frozen, 
calorically perfect thermochemical data to be determined by chemical equili- 
brium calculations. NCEC71 = 2 permits equilibrium data to be determined 
and shifting equilibrium calculations are used for the throat solution. 

NCEC71 = 3 permits frozen flow data to be determined and used for the throat 
solution. The program does not presently allow recalculation of thermochemical 
data during each chamber static pressure iteration loop for MODECH =1 or 2. 
Modification of subprograms REACT and EQULIB would be necessary to reconsider 
reactants data and solution calculation time would be increased considerably. 
When the chamber static pressure is fixed, reactants input data cards used 
for thermochemical data calculations need only be used once by the ther- 
mochemical calculation subroutines since mass flow iterations do not 
affect thermochemical calculations. If thermochemical frozen, calorically 
perfect data is input under namelist EQU, the data remains unchanged for 
iterations for the case. NAORW = 2 before being read by namelist NUM and 
remains unchanged for stacked cases unless read by a namelist NUM read. 

30 
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NAORW = 1 is used when it is desired to specify fixed injector areas and 
NAORW = 2 is used to specify fixed injector velocities. MGORLI = 2 before 
being read by namelist NUM.' If MGORLI = 1 is input, the gas constants 
of gaseous fuel and oxidizer reactants in each zone must be read in an 
initial case by namelist NUM. This option is useful for gaseous injection 
where reactant molecular weights are known. If MGORLI = 2 is input or if 
MGORLI does not appear in an initial case namelist NUM read, fuel and ox- 
idizer reactants densities in each zone must be read in an initial case 
namelist NUM read. MTOTMX = 0 provides that the burned gas in the chamber 
is not to be fully mixed for the throat solution calculations. MTOTMX = 1 
provides that fully mixed flow of the chamber burned gases by used for the 
throat solution calculations. For MTOTMX = 1, one flow zone is used for 
throat solutions and the options NCEC71 = 0 or 1 may be considered. No fully 
mixed solutions may be considered for NCEC71 = 2 or 3. If only one injector 
zone is considered for calculations, MTOTMX = 0 must be used. NSPPG = 0 
before being read by namelist NUM. If NSPPG = 0, special chamber routine 
SPPG is not used. If NSPPG = 1, routine SPPG will be used. NSPPG = 1 may 
only be used for single zone flow of low temperature nonreacting gases. 

For striated flow analysis of rocket motors or for single zone reacting 
gas flows, routine SPPG may not be used. Routine SPPG provides a method 
for cold flow, single zone model analysis found in the laboratory. 

Injector data includes the following: 

(a) Number of flow zones; NCHA1 and NCHA2 

(b) Oxidizer to fuel mass ratios; OXFULE 

(c) Injector reactants temperatures; TIFU and TI0X 
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(d) Injector area data; AITAC, AIPCT, and AIPCTF 

(e) Injector velocities; WIFU and WIOX 

(f) Injected reactants densities; RHOIFU and RHOIOX 

(g) Injected gaseous reactants gas constants; RIFU and RIOX 

Either injector area data or injector velocities and either densities or 
gas constants are input by namelist NUM. Other injector data required 
for initial cases and these also need not be changed for stacked cases 
unless required. 

NCHA1 is required for an initial case read and must be greater than 
zero. NCHA2 = 0 before being read by namelist NUM initially and the sum 
of NCHA1 and NCHA2 cannot exceed twenty. Two regions of flow zones desig- 
nated by NCHA1 and NCHA2 have been used to eliminate or change outer flow 
zones in stacked cases with few changes or to enable the user to modify 
the throat integration technique in routine HPCG for future adaptations. 

Oxidizer to fuel mass ratio is required with all cases and indices 
accounting starts for the inner most flow zone. It may be necessary to avoid 
OXFULE = 0.0 or other very small fractions for a flow zone when using 
chemical equilibrium calculations due to the lower temperature limit on 
combustion products allowed by modified program CEC71 . 

Injector fuel temperature, TIFU, and oxidizer temperature, TIOX, must 
be input in each initial data case with namelist NUM. If OXFULE = 0.0 
for any zone, special handling is not required; however, if oxidizer only 
is injected into a flow zone, the injected material must be designated 
a fuel and OXFULE is correspondingly set equal to zero for that zone. 

If NAORW = 1, AITAC, AIPCT, and AIPCTF are required as initial case 
data by namelist NUM. AITAC is the ratio of total injector areas to the 
chamber area. AIPCT are the ratios of fuel plus oxidizer areas for each 
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flow zone to the total injector area and AIPCTF are the ratios of fuel 
injector area to the fuel plus oxidizer injector area in each flow zone. a 
If AITAC is not input or is zero, AIPCT then must be the injector areas 
for each zone. Thus, injector areas as opposed to area ratios may be 
input directly. 

As an alternative to specified injector areas, fuel injector veloci- 
ties, WIFU, and oxidizer injector velocities, WIOX, may be input. This 
option is used when NAORW = 2. NAORW = 2 before being read by namelist 
NUM for the first case being considered. 

Injector reactant densities, RHOIFU for the fuels, and RH0I0X for 
the oxidizers, are required for MGORLI = 2. RHOIFU (1) and RH0I0X (1) 
correspond to the fuel and oxidizer densities in zone one, etc. MGORLI = 2 
before the first case data are read by namelist NUM. 

As an alternative to injector densities, fuel gas constants, RIFU, 
and oxidizer gas constants, RIOX, may be input for gaseous injection. 

The option is exercised by setting MGORLI = 1 in namelist NUM. 

Refer to Appendix D for more information on data read by namelist 
NUM and to Appendix F for examples. 


a If MODECH = 2 an exception to injector area input must be made. 
AITAC must be the ratio of fuel injector areas to the total chamber area 
and AIPCT is the ratio of fuel injector area for the ith zone to the 
total fuel injector area. For MODECH = 2, AIPCTF is not required input 
but is internally set equal to 1.0. 
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Namelist EQU 


If NCEC71 = 0, frozen, calorically perfect thermochemical data is used 
for chamber and throat calculations. The data required for each flow zone 
are the total temperature, TOG, the isentropic coefficient, GAMG, the burned 
gas molecular weight, MWG, and an integer variable NCHANG. Each case 
requires one set of data read by namelist EQU when NCEC71 = 0 and the data 
set follows variables read by namelist NUM. Each initial case requires 
all thermochemical variables described above and NCHANG must be greater 
than zero. For stacked cases where thermochemical data does not change 
from the previous case, it is necessary only to input NCHANG = 0. Table I 
is a schematic of the input and namelist variables required for all problems 
when NCEC71 = 0. Operation of the program with this input consititutes 
the first kind of option mentioned at the beginning of this section. 

THERMO Data 


When thermochemical data are to be determined by chemical equilibrium 
calculations, NCEC71 = 1, 2, or 3, two types of data are required. The 
first type of chemical data required for this kind of problem as indicated 
by Table II, are THERMO Data. This consists of a THERMO code card, tempera- 
ture range data, and a library of coefficients for calculating thermodynamic 
parameters of chemical species. The data used by this program is identical 
to that used by NASA SP-273 . ^ The data may be read by card, tape, or disk 
and is written on unit 4 for use during a run. A permanent tape or disk 
containing the data may be made during a run (requiring a THERMO code card) 
or prior to a run (without THERMO code card) by proper job control cards. 
Options and changes to the THERMO Data of NASA SP-273^ are possible with 
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this program also. The authors have found it convenient to permanently 
store THERMO data on a disk and use the disk with this program or with 
NASA program CEC71 . 

REACTANTS Data 

If NCEC71 = 1, 2, or 3, complete sets of REACTANTS Data must be input 
for each initial case. There must be as many REACTANTS data sets as there 
are flow zones; however, for a stacked case after an initial case where 
the thermochemical data will not change, the REACTANTS data sets need not 
be input. If cases are to be stacked as indicated above, the Control Card, 
LETOUT, must have the integer 71 starting in column 4 for each stacked 
case. The REACTANTS data set consists of three types of cards. The first 
type in this set has REACTANTS punched in columns one through nine. The 
second type of cards are reactant cards coded with the same format as the 
reactant card coding found in NASA SP-273^ except that only 0/F ratios may 
be used. A blank card follows the reactants cards. The third type of card 
in each REACTANTS set has NAMELISTS punched in columns one through nine. 
This card initializes chemical equilibrium calculations. 

Additional information on program inputs and options may be seen in 
Appendixes D and F. 

Units 

In this program the English Technical system of units is used with 
the exception of thermo data and the reactant card inputs which are used 
by the chemical equilibrium calculations. Thermo data and reactant card 
inputs form an optional kind of input and do not use the English Technical 


35 



system so as to be compatible with the inputs currently used by NASA program 
CEC71. The English Technical system is as follows: 


Physical Quantity Unit Symbol 

Length feet ft 

Time second sec 

Temperature Rankine R° 

Force pound Ibf 

Pressure pounds per psf 


square foot 

Mass Flow pound seconds lbf-sec/ft 

per foot 

The Universal gas constant used is 1545.43 Ibf- ft/1 bmole-°R. Density 

2 4 

has the units of lbf-sec /ft and specific impulse has the units of 

lbf-sec/lbmass. Where required, a conversion constant of 32.174 lbmass-ft/ 

2 

lbf-sec is used. The oxidant to fuel weight ratio is used exclusively if 
reactant cards are used. Symbols used to represent parameters may be seen in 
Appendix B and program variable names are given in Appendix C. 
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DESCRIPTION OF PROGRAM OUTPUT 


The final output of the program is in the form of printed tables 
that are designed to be self-explanatory. Except where fully mixed flow 
is considered all output are the same. 

Namelist NUM 

For an initial case, parameters in namelist NUM may be printed by 
a namelist write. The printout option is exercised by setting LETOUT = 1. 
For LETOUT >_ 2 namelist NUM variables are not printed. 

Heading and Comments 

The case number is printed for each case. The numbers are inter- 
nally generated and not user controlled. Information on the input comment 
card are printed following the case number. 

Option Data 

Option data includes MODECH, MGORLI , NAORW, NCEC71 , NSPPG, MTOTMX, 
LETOUT, and NCHA1 + NCHA2 or the number of zones designated by NINCH. 

Also included are the number of throat solution iteration loops required 
for convergence (LOOP) and R c /R.j.. 

Geometry and Pressure Output 

This output includes R^, the converged value of the wall radius at 
the throat as determined by the throat solution, and the chamber contraction 
ratio, EPSC. R T in this output is the specified throat radius. Pressure 
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data includes the chamber static pressure, P , and the centerline static 
pressure at the throat, P^. 

Injector Parameters 

Output includes the zone stream function limits, the flow zone index, 

Pj » Tt , T| , Aj , Aj , Wr , Wy , Py , P y , Rt , Ry , 

fu^ ox^ fu^ ox^ fu.j ox^ rUy- ox^. *fu.j 1 ox i 

• • 

m T , m T , P R and P R 
A fu. 1 ox i K f Ui ox^- 

If MODECH = 3, reservior pressures are set equal to zero and if MGQRLI = 2 
injector gas constants are set equal to zero. If 0/F i is zero, oxidizer 
parameters are zero. 

Burned Gas Parameters 

When MTOTMX = 0, striated flow at the thraot is considered and com- 
bustion product output parameters for the chamber included P , T , 

% °9i 

Cp » Y a » T » A , w , p , and m. for each flow zone, 
g. y i y i 9 i y i 9 i 1 

Mixed Gas Parameters 

For fully mixed flow in the chamber (MTOTMX > 0) only one set of burned 
gas parameters are output. These parameters are designated by the same 
symbols as were zone parameters in burned gas outputs. In addition, however, 
the chamber static pressure of the fully mixed gases and the total pressure 
of the burned gases in each flow zone prior to mixing are included. 
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Throat Plane Output 


Spatial distributions of r, P, ij>, T, w, M, and y are included 
in throat plane output. Values at the throat wall for the converged 
solution are also given in the output. 

Performance Output 

Vacuum specific impulse at the throat, characteristic velocity, 
total mass flow rate for radius of curvature effects and one 
dimensional flow are performance outputs. In addition, the ratios of mass 
flows (discharge coefficient), vacuum specific impulses, and characteristic 
velocities are printed. Output units are given. 

Units for all outputs which appear in the other outputs are found in 
the UNITS section of this report. 
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PROGRAM MODULES 


The program consists of two modules as depicted in Appendix G. The 
first module was developed specifically for striated nozzle flow with 
throat radius of curvature effects. The second module is an adaption 
of NASA Lewis program CEC71 documented in NASA SP-273^ . The original 
version of CEC71 has been modified to be compatable with requirements 
of the first module. 

Striated Nozzle Flow 

This module consists of subprograms which read comments, control data, 
and namelist data. If subprogram DRIVCE is made a dummy subprogram, the 
striated nozzle flow module may be operated without chemical equilibrium 
subprograms for frozen, calorically perfect gas flows. In this form, the 
program would consist of a main program and fourteen subprograms: two 

of which, L0AD1 and TAB, may also be made dummy subprograms. The purpose 
of this module is to determine and write injector, chamber burned gas, and 
throat parameters. 

Chemical Equilibrium Calculations 

This module, depicted in Appendix G, is used to determine thermochem- 
ical parameters for use by the first module. Linkage between the first 
module and this module is by common block and conversions within subprogram 
LOAD! and entry to L0AD1 , L0AD2. This module consists of sixteen subprograms 
and a block data subprogram. General output of program CEC71 has been 
suppressed; however, error message output has been retained. 
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The chemical equilibrium calculations module reads THERMO data and REACTANTS 
data inputs with formats and units of the unmodified form of program CEC71 . 
Removal of comments from data output cards in modified CEC71 in this program 
would cause data output to appear after the heading and comments unless 
a stacked case with LETOUT =71 is being run. 
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ROUTINES 


The first module of this program will be described in this section 
whereas Reference 1 describes in detail routines of module two (modified 
CEC71). Limitations on inputs to module two are described in the DESCRIP- 
TION OF PROGRAM INPUT section. 


Main Program 

A flow chart of the main program is given in Appendix H. Generally, the 
routine performs the following functions: 

(1) Initializes inputs and constants. 

(2) Reads namelist NUM. Writes NUM if LETOUT = 1. 

(3) Calculates chamber geometry. 

(4) Initializes P g for MODECH = 1 or 2. 

One dimensional equations are used to estimate initial values of P g 


or m.j . For fixed mass flow in zones, initial estimates of , 

MWg , and T q are required. For fixed reservior pressures, the cham- 
ber pressure is initially estimated to be a fraction of the fuel 
reservior pressure in the first flow zone. The chamber static pres- 
sure is initially estimated if required in the main program. 

(5) Calls subroutine THROAT. 

The purpose of the main program is to serve as a driver for the remain- 
der of the program. 

Subroutine THROAT 

The purpose of subroutine THROAT is to control the throat plane itera- 
tion technqiue, determine performance parameters, and write results. A flow 
chart of subroutine THROAT is given in Appendix H. This subroutine performs 
throat plane calculations prior to performance parameter calculations and so 
the two types of calculations may be discussed separately. 
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Initially, the number of iterations performed to satisfy choking at the 
throat is checked to determine whether more than seven iterations have been 
performed to determine the total mass flow or chamber static pressure re- 
quired. If LOOP = 8 then results of the seventh loop are used to determine 
output parameters. If mass flows are specified or if reservoir pressures 
are specified, P g is being sought in the throat solution. For these modes, 
MODECH = 1 or 2, alterations of P g on the seventh loop is made by interpola- 
tion when Pg has been bracketed by previous iteration values of the chamber 
static pressure. If it has not been bracketed, the sixth iteration is used 
for output parameters. If specified reservoir pressures are considered, an 
additional requirement is that Pg must be less than one percent smaller than 
any fuel reservoir pressure. Very seldom more than four iteration loops are 
required for a converged throat solution. 

The chamber solution is determined and then throat plane centerline 
static pressure estimated and altered to determine the throat radius corres- 
ponding to maximum mass flux. Integrations performed by subroutine HPCG 
require boundary conditions for both dependent variables at the nozzle center- 
line. These are P = P^ and 4) = PSIOl(l) = 0.0. After itegration to 4>=1.0 
is complete, the corresponding wall radius is determined by subroutine AITKEN. 
Radii are ordered monatonically upward by subroutine UPORD. Successful deter- 


mination of a minimum radius consists of two criteria which permit checks to 


be made as to (1) whether the radius is a minimum based on (dr/r) < .05% and 

w 

dP CL^ P CL <1 ^ and ^ whether the minimum radius determined is within .05% of Rj. 


If the first of the above is satisfied but not the second. 


P 9 


or m 


total 


is altered and iteration starts again by indexing LOOP, Subroutine UPORD is 


used to order the P CL during a loop once the minimum radius is bracketed for 
that loop. A second order curve fit by subroutine PLYCF is made to determine 
P CL for ^ dr w /,dP CL^ = 0,0 ‘ The first criteria stated above must be satisfied 
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or the process is repeated with the three smallest radii of the loop. 

With the first and second criteria stated above satisfied, remaining 
parameters except throat plane parameters and performance are output. 
Performance calculations are made as described in Performance Solution. 
Throat plane and performance parameters are printed as described in Throat 
Plane Output and Performance Output, respectively. 

Subroutine CHAMBR 

The purpose of subroutine CHAMBR is to determine chamber parameters 
for striated or fully mixed flow. A flow chart for subroutine CHAMBR is 
given in Appendix H. This routine performs calculations shown in Chamber 
Solutions and calls other subroutines which cause thermochemical data to 
be calculated or read. 

Upon being called by subroutine THROAT, subroutine CHAMBR initializes 
the iteration indexes NOAI and LOOPC and begins iteration loops by initi- 
alizing the injector static pressure, Pj. For MODECH = 2, subroutine PRSPFD 
is called to determine zone mass flows and oxidizer injector areas. Both 
fuel and oxidizer reservoir pressures are input and zone 0/F ratios must 
be specified. If MODECH = 3, subroutine P6FIX is called to determine 
thermochemical data and to determine zone mass flows from the total mass 

flow. Chamber parameters T , A , and w are next intialized for the 

9^ =i 9i 

iterations and, for MODECH = 1 or 2, thermochemical data are determined 
by a direct call to subprogram EQULIB. Special routine SPPG is called 
to provide the chamber parameters if option parameter NSPPG = 1 is in 
effect. Injector parameters are determined and chamber burned gas para- 
meters T , p , w , and A for iteration loop LOOPC are calculated. 

9i 9i 9i °i 

The primary iteration with index NOAI is used for M0DECH=2 which, after a 
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complete iteration with index LOOPC, alters the mass flow by repeated calls 
to PRSPFD. Index NOAI also serves to control the injector pressure iteration 
calculations and a convergence error message is provided when NOAI > 20. For 
MODECH = 2, iteration is complicated by the fact that mass flow must satisfy 
reservoir equations and chamber equations. Thus, for specified reservoir 
pressures, one equation set is added and one set of unknowns, m^ , must exist 
to determine a chamber solution. 

Convergence of the chamber solution permits dividing streamlines to be 
determined and corresponding radii in the chamber and total pressures to be 
determined. If MTOTMX > 0, fully mixed flow parameters are determined by 
equations given in the section Fully Mixed Flow. 

Subroutine PGFIX 

This subroutine is called by CHAMBR when MODECH =3. The purpose of 
this routine is to call subrotine EQULIB for thermochemical parameters. It 
also serves to provide an initial estimate of the total and zone mass flows 
and to calculate zone mass flows during iterations with index LOOP. 

Subroutine PRSPFD 

This routine provides an initial estimate for Pj for NOAI =1. It also 

provides intermediate calculations for zone mass flows and oxidizer injector 
areas. PRSPFD is used when MODECH = 2. 

Subroutine SPPG 

Subroutine SPPG is used for single zone low temperature or non-reacting 
gas flow chamber calculations when MODECH = 1 or 3. Its use is optional and 
occurs when NSPPG = 1. The iteration index LOOPSP permits simultaneous solution 
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of chamber equations to satisfy convergence requirements. Convergence 
requirements are that burned gas static temperature for loop LOOPSP + 1 
be within .05% of the temperature of loop LOOPSP. 

Subroutine EQULIB 

This routine serves three purposes. These are: 

(1) Read frozen, calorically perfect thermochemical data under name- 
list NUM when NCEC71 = 0. 

(2) Call subprogram DRIVCE (formerly the main program for CEC71 ) for 
NCEC71 = 1, 2, or 3. 

(3) Store chamber thermochemical data for reuse during iterations. 
This is necessary when MTOTMX = 1. 

Logical IF statements in subprogram EQULIB determine the order of 
operations so to satisfy option and iteration loop requirements. 

Subroutines L0AD1 and L0AD2 

Subroutine LOAD! (L0AD2 is an ENTRY to LOAD!) contain COMMON blocks 
POINTS and PERF which permit data transfer. These two blocks appear 
in NASA Lewis program CEC71 . In addition, COMMON blocks ZONE, ONE, and 
ACSTAR have been included to permit data transferral between the two cal- 
culation modules which make up this program. L0AD1 and L0AD2 apply con- 
version units and store thermodynamic parameters for chamber and throat 
calculations. Each flow zone is assigned on index, KZONE, and two di- 
mensional arrays permit interpolation in a flow zone at the throat for an 
internally generated pressure schedule. Ten pressure ratios are generated 


46 



for each flow zone in subprogram ROCKET when NCEC71 = 2 or 3. In addition, 
one dimensional sonic pressure and chamber pressure are stored. Thus, 
thermochemical calculations are performed for twelve pressures for each 
flow zone when NCEC71 = 2 or 3. When NCEC71 = 1, only data corresponding 
to the chamber and sonic pressure are determined; however, data corres- 
ponding to sonic pressure are not used for calculations ince NCEC71 = 1 is 
for frozen, calorically perfect flow. 

Subroutine HPCG 

This subroutine uses Hanning's modified predictor-corrector method 
for the solution of throat plane equations which are of the initial -value 
type. This subroutine is discussed and listed with operating instructions 
in Reference 9. HPCG uses two external subprograms, FCT and OUTP. FCT 
provides the right hand side of two first order equations and OUTP stores 
the solution vectors of HPCG. Few changes have been made to the original 
version of HPCG. Comment cards have been removed and statements which 
have been added are so indicated by alphabetic characters in card sequences 
replacing last character zeros. Look to the program listing of HPCG to 
see the changes. 

Subroutine FCT 

Subroutine FCT provides the right hand side of the two first order 
ordinary differential equations for P and with r as the independent vari- 
able. In the call to FCT from HPCG three parameters are transferred; R 
is the independent variable, Y are the solutions P and ip, and DERY are the 
first two error weights for the dependent variables which are saved in 
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THROAT under variable name ERWT. FCT provides a check of option NCEC71 
to determine if frozen, cal ori call y perfect flow is being considered or if 
frozen or equilibrium flow is being considered. Determination is made as 
to which flow zone integrations are being performed in and provides the 
corresponding thermochemical variables. For NCEC71 = 0 or 1 , the index K 
determines the thermochemical variables to be used. For NCEC71 = 2 or 3, 
subroutine TAB is called to interpolate on the current value of K and Y(l). 
Y(l) is the current pressure, P, and K corresponds to KZONE. KZONE is 
generated in DRIVCE and is the index on two dimensional arrays developed 
in L0AD1 . Equations used for dependent variables are described in the 
THEORY section under Throat Solution. 

Subroutine OUTP 

Subroutine OUTP is called by HPCG. This routine stores dependent 
variables with index J. If NCEC71 = 0 or 1 , temperature (T ( J ) ) Mach number 
(YPSAV(J,4)), and velocity ( YSAVE( J ,3) ) are determined in OUTP. OUTP 
also causes PRMT (5 ) to be set equal to 1.0 when ij> ( YSAVE( J ,2) ) becomes 
>_ 1.0. PRMT(5) f 0.0 causes return to subroutine THROAT from subroutine 
HPCG. 

Subroutine TAB 

TAB is an interpolation routine used by subroutine FCT for deter- 
mination of thermochemical data when NCEC71 = 2 or 3. 
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Subroutine AITKEN 


This subroutine is called by subroutine THROAT to determine the 
throat wall radius after return to THROAT from HPCG. It also provides 
wall values of throat plane parameters which are determined in the performance 
section of subroutine THROAT. 

Subroutine UPORD 

This subroutine is called by subroutine THROAT to order a pair of 
real array variables. The parameter list for this subroutine contains 
two array variables and two integers. The first integer gives the number 
of pairs and the second integer, MORDER, determines which array will be 
ordered up. MORDER f 0 orders the first array upward with corresponding 
variables in the second carried along. MORDER = 0 orders the second array. 
MORDER = 1 for this program. 

Subroutine PLYCF 

Subroutine PLYCF is a least squares polynomial curve fit routine called 
by subroutine THROAT to predict p~. for r 

w min 

Subroutine DRIVCE is the modified main program of NASA Lewis program CEC71 . 
Subroutines REACT, SEARCH, HCALC, SAVE, EQLBRM, CPHS, MATRIX, GAUSS, 0UT1 , 
VARFMT, EFMT, THERMP, ROCKET, RKTOUT, FROZEN, and BLOCK DATA are subprograms 
used by the program CEC71 and are described in NASA SP=273^ . Differences 
between these routines and the routines of CEC71 may be seen by card se- 
quence numbers. A / symbol after the sequence name indicates a comment 
card has been made and an alphabetic character and a repeated sequence 
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number after the sequence name indicates a statement has been added. 

Changes to CEC71 occur in DRIVCE (CEC71 program Main), ROCKET, and RKTOUT. 
Subroutines SHCK and DETON of CEC71 have been removed for this modification. 
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APPENDIX A 


PHILOSOPHY 

The philosophy used in developing this program has been to define 
where possible, subprograms or variables with names which are appropriate 
and representative of the purpose they serve or the commonly used symbol 
they represent, respectively. Use of the same variable name to represent 
different variables has been avoided. Subprograms which have been written 
and documented by others which serve as utility subprograms may have changes; 
however, changes which have been necessary have been made so as to be easily 
detectable. Individual card sequence numbers expose the changes made. 

The program has been written to be as machine independent as possible. 

The FORTRAN IV language (FORTRAN G level 20) has been used and the program 
is constructed so that user modifications are possible, particularly where 
a change in the units for input and output is desired. Transferral of var- 
iables between subprograms is by common block where possible and algorithms 
have been used which are logical and easy to follow. Computation time for 
typical cases are small (approximately one minute on the IBM 360/65). The 
main input data set is by namelist to simplify checking input cards and 
to simplify stacked case input data sets. 

The program is intended to provide a capability for rapid analysis 
with detailed chamber operating conditions and spatial distributions of 
thermodynamic flow properties on the nozzle throat plane. An important 
aspect of the program is that it permits throat plane caluclations with a 
variable adiabatic or isentropic coefficient and is not limited to a constant 
stagnation pressure in the chamber. The computer program for calculation 
of chemical equil ibrium compositions as found in NASA SP-273^ has been 
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incorporated as subprograms. Changes to this program, CEC71 , may easily 
be seen in the program listing and compaired to the original version of 
the program by the sequence code found in the last eight columns of the com- 
puter card. The SHCK and DETON subprograms of CEC71 were removed and re- 
ferences to them have been made comment cards. Input data to the modified 
form of CEC71 in this program are simplified somewhat but, where required, 
have the same form and units used by program CEC71 . Standard output of 
the modified CEC71 calculations have been suppressed by commenting output 
in the manner described above; however, error messages from this program 
have not been suppressed. 
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APPENDIX B 


A 

A. 


C 

C r 


9 C 

H o 

h 


sp 


vac 


M 


MW 


m 

0/F 


P 


R 


R, 

<3, 


r 

T 


SYMBOLS 

2 

area, ft 

2 

chamber area, ft 
2 

throat area, ft 

P/tp) Y , poly tropic equation of state constant (isentropic) 
discharge coefficient 

2 2 

constant pressure specific heat, ft /sec -°R 

characteristic velocity, ft/sec 

32.174 lbmass-ft/lbf-sec 2 

2 2 

total enthalpy, ft /sec 

2 2 

enthalpy, ft /sec 

throat vacuum specific impulse, lbf-s.ec/lbmass 
Mach number 

molecular weight, Ibmass/lbmole 
mass flow rate, Ibf-sec/ft 
oxidant to fuel mass ratio 

2 

stagnation pressure, lbf/ft 
2 

pressure, lbf/ft 

2 2 

gas constant, ft /sec -°R 
throat wall radius of curvature, ft 
throat radius, ft. 

flow radius of curvature at throat, ft 
universal gas constant, 1545.43 lbf-ft/lbmole-°R 
radius; radial coordinate, ft 
temperature, °R 
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stagnation temperature, °R 
radial velocity component, ft/sec 
tangential velocity component, ft/sec 
axial velocity component, ft/sec 
axial coordinate, ft 

isentropic exponent or adiabatic coefficient 


3.14159265 

density, lbf-sec^/ft^ 
stream function = 


m(r) 

/ ['“'"'total] 


centerline value on or near the throat plane 


combustion product or burned gas 

injector 

minimum value 

fully mixed 

oxidizer 


reservoir 


total or whole amount 


wall value on or near the throat plane 
one dimensional 


Superscripts: 


modi fied 

critical throat condition 


Indices: 


striation zone (inner most zone is 1) 
local value on the throat plane 
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APPENDIX C 


COMMON VARIABLES 


All common variables are given in this appendix except those developed 
in the original unmodified program CEC71 ,which is documented in NASA SP-273. 


Variable 

Dimension 

Common 

Label 

Description and Comments 

AC 

1 

THR 

A c • 

AG 

20 

THR 

A g * 

AIFU 

20 

THR 

A l * 

fu i 

AIOX 

20 

THR 

A I * 

ox i 

AIPCTF 

20 

THR 

Fraction of injector fuel 
area assigned to each zone. 

AMAC 

20,13 

LOTAB 

Mach number array in each zone 
corresponding to APRES values. 

APRES 

13 

LOTAB 

Pg/P array used for interpolation 

at the throat plane for frozen 
or equilibrium flow. 

ARHO 

20,13 

LOTAB 

Density array in each zone 
corresponding to APRES values. 

ASON 

20,13 

LOTAB 

Sonic velocity array in each 
zone corresponding to APRES values 

AUX 

16,4 

blank 

Auxiliary storage array (See HPCG 
subprogram description). 

AVEL 

20,13 

LOTAB 

Velocity array in each zone cor- 
responding to APRES values. 

C 

1 

CALC 

P g^ p g f 9 i ’ used for f rozen ’ ca ^ _ 

orically perfect flow only. 

CDINJ 

20 

THR 

Injector discharge coefficients 
for MODECH = 2. 


a For MODECH = 2, AIPCT becomes the ratio of fuel injector areas to chamber 
area and AIPCTF is not required. 
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CLP 

10 

bl ank 

CPG 

20 

TWO 

CST 

20 

LOAD 

CSTAR 

13 

ACSTAR 

DERY 

4 

blank 

DIACH 

1 

THR 

EPSC 

1 

ONE 

ERWT 

4 

blank 

GAM 

1 

CALC 

GAMG 

20 

TWO 

GAMR 

20,13 

LOTAB 

gamsav 

200 

CALC 

HO 

(Hzero) 

1 

CALC 

ICASE 

1 

ONE 

IHLF 

1 

ONE 

J 

1 

TWO 

K 

1 

CALC 

KZONE 

1 

ZONE 



c*.j from chemical equilibrium 
calculations. 

c* corresponding to P g for each zone 

Input error weights, dependent 
variable derivatives (see FCT and 
HPCG subprograms descriptions). 

Chamber diameter. 

Chamber contraction ratio. 

Error weights for dependent vari- 
ables in subprogram HPCG. 

Local value of y on the throat 
plane. 



Isentropic coefficient array in 
each zone corresponding APRES values 

y array on the throat plane. 

Local value of total enthalpy, H q 



Solution case index. 

Step size bisection parameters used 
by subroutine HPCG. 

Step index for local throat plane 
parameters 

Index parameter for zone designa- 
tion . 

Zone index. 
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LETOUT 

1 

ZONE 

Integer control constant. See 
input data for LETOUT options. 

LOOP 

I 

blank 

Index counter designating 
the number of throat solution 




iteration loops. 

1 

MDTI 

20 

THR 


MDTIFU 

20 

THR 

« 

m T 

A fu. 

i 

MDTI OX 

20 

THR 

rrij 

ox. 

MGORLI 

1 

ONE 

Input option parameter indicating 
injector gas constant input (=1) 




or injector density input (=2). 
MGORLI = 2 initially. 

MODECH 

1 

blank 

Input option parameter indicating 
the mass injector scheme. See 
input data for MODECH options. 

MORDER 

1 

blank 

= 1, describes the method for 
ordering variables in subroutine 
UPORD. 

MTOTMX 

1 

THR 

Input option parameter. MTOTMX 
= 0 for striated flow, MT0TMX=1 
for fully mixed flow. MT0TMX=0 
initially. 

MWG 

20 

THR 

MW . 

9i 

MO 

(Mzero) 

1 

TWO 

2m., total mass flow. 

• 

MOCH 

MzeroCH) 

7 

blank 

2m.. corresponding to LOOP. 

NAORW 

1 

THR 

Input option parameter indicating 
whether injector areas are fixed (=1) 




or injector velocities are fixed 
(=2). Initially NAORW = 2. 

NCEC71 

1 

ONE 

Input option indicating the degree 
of chemical equilibrium. See 
Input data for NCEC71 options. 
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NCHA1 

1 

TWO 

The number of flow zones in an 
arbitrary chamber region (1<NCHA1<10 
NCHA1 is an input option. 

NCHA2 

1 

TWO 

The number of flow zones in an 
arbitrary chamber region (0<NCHA2<10 
NCHA2 is an input option. 

NDIM 

1 

ONE 

= 2. See HPCG subprogram description 

NSPPG 

1 

ONE 

Special input option for single 
zone nonreacting flow calculations. 
Initially NSPPG = 0. See input list. 

OXFUL 

20 

ZONE 

0/F.j for chemical equilibrium cal- 
culations . 

OXFULE 

20 

THR 

0/F. . 

PCTMO 

(PCTMzero) 

20 

THR 

Fraction of total flow in the ith^ 
flow zone. 

PG 

1 

THR 

V 

PGCH 

7 

blank 

Pg corresponding to LOOP. 

PGPSIA 

1 

ZONE 

Pg internally converted to psia for 
chemical equilibrium calculations. 

PI 

20 

THR 

p i 

PMIX 

1 

MIXED 

P 9 ■ 
a mix 

PRF 

20 

THR 

P R f ‘ 
fu i 

PRO 

20 

THR 

P R • 

ox-i 

PRMT 

6 

bl ank 

Input and output array for HPCG 
(see input desription). 

PSIOl 

10 

TWO 

*. in the first flow region. 

PSI02 

10 

TWO 

a. in the second flow region. 
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PS112 

1 

TWO 

ij> corresponding to the streamline 
between the first and second flow 
regions . 

PO 

(Pzero) 

20 

MIXED 

Total pressure in each zone before 
full mixing occurs. 

POOL 

(PzeroCL) 

7 

blank 

P Q corresponding to LOOP. 
9 i 

POG 

(PzeroG) 

20 

TWO 

p o • 

RAD 

20 

THR 

Radius corresponding in the chamber. 

RCRT • 

1 

ONE 

R c /R t . 

RCWALL 

1 

TWO 

R c • 

RGAS 

20 

THR 

\ • 

RHOG 

20 

THR 


RHOIFU 

20 

THR 

P I 

fu i • 

RHOIOX 

20 

THR 

Pj 

ox i 

RIFU 

20 

THR 

R 

fu i 

RIOX 

20 

THR 

R I 

OXi 

RUM 

1 

blank 

not used. 

RMIN 

7 

bl ank 

r w corresponding to LOOP. 

RMINSQ 

7 

bl ank 

2 

r corresponding to LOOP. 

RT 

1 

TWO . 

r t • 

RW 

10 . 

bl ank 

Wall radius vector corresponding 


to LOOP. 
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T 

200 

ONE 

T on the throat plane. 

TEMPE 

20,13 

LOTAB 

Temperature array in each zone 
corresponding to APRES values. 

TG 

20 

THR 

\ * 

TIFU 

20 

THR 

T I 

fu i 

TI0X 

20 

THR 

T I * . 

ox i 

TWO PI 

1 

TWO 

2ir 

TOG 

(TzeroG) 

20 

THR 

T o • 

g. 

UNIVR 

1 

THR 


URMNSQ 

7 

blank 

? 

1/r corresponding to LOOP. 

VAC1DI 

20 

LOAD 

>'sp • 

K vac 

VACSPI 

20,13 

LOTAB 

Vacuum specific impusle array 
in each zone corresponding to 
APRES values. 

WG 

20 

THR 

\ ' 

WIFU 

20 

THR 

W I 

fu i 

WIOX 

20 

THR 

W I 

ox-i 

WMR 

20,13 

LOTAB 

Molecular weight array in each 
zone corresponding to APRES values. 

XSAVE 

200 

TWO 

r on the throat plane. The index 
is J. 

Y 

4 

bl ank 

Y(1)=P and Y(2)=ip on the throat 


plane corresponding to XSAVE. 
Y(3) and Y(4) are not used. 
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YPSAV 

200,4 

TWO 

YPSAV(J,1) is the pressure derivative 
array on the throat plane. YPSAV (J,2) 




is the derivative of ^ , 

YPSAV ( J ,3 ) is T q /T on the throat 

9i 

plane. YPSAV(J,4) is the Mach 
number array on the throat plane. 

YSAVE 

200,4 

TWO 

YSAVE (J ,1 ) = p on the throat 
plane. YSAVE(J,2) = ij> on the 
throat plane. YSAVE(J,3) = w on 
the throat plane. YSAVE(J,4) 
is dA corresponding to YSAVE (J , 1 ) - 

ZCH 

1 

THR 

not used. 
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DIMENSIONED VARIABLES 


Variable 

Dimension 

Routines 
where used 

Description and Comments 

A 

11 

PLYCF 

polynomial coefficients. 

ACO 

11 

THROAT 

polynomial coefficients. 

AGAMGU 

20 

MAIN 

initial estimates. 

AIPCT 

20 

MAIN 

fraction of total injector 
area for zone i . a 

AMACH 

200 

THROAT 

Mach number on the throat plane 

AMWGU 

20 

MAIN 

initial MW estimates. 

ATOGU 

(ATzeroGU) 

20 

MAIN 

initial T estimates. 

AUX 

16,2 

HPCG 

Auxiliary storage array (see 
HPCG subprogram description). 

CLP 

10 

UPORD 

p 

r CL* 

DERY 

2 

FCT, OUTP, 
HPCG 

input error weights, dependent 
variable derivatives (see FCT 
subprogram description). 

DRDPSV 

10 

THROAT 

dr w /dP a . 

DUM1 

3 

THROAT 

independent dummy variable for 
AITKEN interpolation. 

DUM2 

3 

THROAT 

dependent dummy variable for 
AITKEN interpolation. 

PII 

20 

CHAMBR 

stored P. array during CHAMBR i 
ation. 1 

PRMT 

6 

OUTP, HPCG 

input and output array for 
HPCG (See input description). 

PS I 

200 

THROAT 

\|i on throat plane. 


a For MODECH = 2, AIPCT is the ratio of fuel injector area to total fuel 
injector area for the ith zone. 
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RADT 

20 

THROAT 

radius of stream tube shell 
boundary at the throat plane. 

RW 

10 

UPORD 

V 

SAVGAM 

20 

EQULIB 

y , values stored to prevent 
9 i 




destruction in fully mixed flow 
cases. 

SAVMWG 

20 

EQULIB 

MW values stored to prevent 
9i 




destruction in mixed flow cases. 

SAVTOG 

(SAVTzeroG) 

20 

EQULIB 

T q values stored to prevent 
9 i 




destruction in fully mixed flow 
cases . 

SUMAGG 

20 

CHAMBR 

2A 

9i 

TGA 

15 

SPPG 

T 9l 

TGA 

20,10 

CHAMBR 

T »i ' 

WHAT1 

20 

MAIN 

first input comment card required 
for each case (20A4) . 

WHAT2 

20 

MAIN 

second input comment card required 
for each case (20A4). 

X 

Y 

variabl e 
2 

AITKEN, 

PLYCF 

FCT,OUTP, 

HPCG 

general independent varialbe. 

dependent variable in HPCG (see 
FCT subroutine description). 

Y 

variable 

AITKEN, 

PLYCF 

general dependent variable. 
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APPENDIX D 


INPUTS AND OPTIONS 


Comment Cards (2): 


Order Contents Format 

First 9 Any information 20A4 

pertinent to the 
case being run. 

Control Card (Variable name LETOUT) ; 

Order Contents Format 

Second* 5 Integer Control 15 

(LETOUT in Common 
label ZONE) 


Card Columns 
1 to 80 


Card Columns 


1 to 5 


(1) If LETOUT = 0, program computations stop and program output begins. 

(2) If LETOUT = 1, any initial case input follows and variables in name- 
list NUM will be output by namelist write. 

p d 

(3) If LETOUT > 1, but f 71, any initial case or a stacked case using 
input thermochemical data follows. Variables in namelist NUM will 
not be output. 

(4) Stacked cases which require thermochemical calculation results be used 
for solutions as opposed to input thermochemical data must have 
LETOUT =71 in columns 4 to 5. This stacked case option is used as 
when the chamber static pressure, flow zone 0/F ratios, and reactants 
data for each flow zone will not change for the stacked cases. Nozzle 
geometry and/or propellant distributions may be changed when LET0UT=71 
for stacked cases which meet requirements stated above. 


a Two comment cards are required for all cases. 

^The control card must be preceded by two comment cards. 

c More than one initial case may be run at one time. Any initial case is a 
case where a complete set of thermochemical input data is required. 

d Stacked cases refer to cases where thermochemical data need not be cal- 
culated (NCEC71 = 1, 2, or 3) or read (NCEC/l = 0). 
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Namelist NUM Variables: 


Variable 

Dimension 

Type 

Common Label 

Value before 

PRMT 

6 

R 

blank 

by first case 
undefined 

ERWT 

4 

R 

blank 

0 . 

PG 

1 

R 

THR 

0 . 

MDTI 

20 

R 

THR 

0 . 

RT 

1 

R 

TWO 

undefined 

RCRT 

1 

R 

ONE 

undefined 

EPSC 

1 

R 

ONE 

undefined 

OXFULE 

20 

R 

THR 

0 . 

PRF ;PRO 

20 ;20 

R 

THR 

0 . 

ESPG 

1 

R 

- 

.7 

AITAC 

1 

R 

- 

0 . 

AIPCT 

20 

R 

- 

0 . 

AIPCTF 

20 

R 

THR 

0 . 

RHOIFU 

20 

R 

THR 

0 . 

RIFU 

20 

R 

THR 

0 . 

TIFU 

20 

R 

THR 

0 . 

WIFU 

20 

R 

THR 

0 . 

RHOIOX 

20 

R 

THR 

0 . 

RIOX 

20 

R 

THR 

0 . 

TIOX 

20 

R 

THR 

0. 

WIOX 

20 

R 

THR 

0 . 

PCTMO 

20 

R 

THR 

0 . 

AGAMGU 

20 

R 

- 

0 . 

AMWGU 

20 

R 

- 

0 . 

ATOGU 

20 

R 

- 

0 . 

CDINJ 

20 

R 

THR 

.7 
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Namelist NUM Variables (cont'd.) 


Variable 

Dimension Type 

Common label 

Value before read by 






first case 

MTOTMX 


1 

I 

THR 

0 

MODECH 


1 

I 

blank 

undefined 

NAORW 


1 

I 

THR 

2 

MGORLI 


1 

I 

ONE 

2 

NCHA1 


1 

I 

TWO 

undefined 

NCHA2 


1 

I 

TWO 

2 

NSPPG 


1 

I 

ONE 

0 

NCEC71 


1 

I 

ONE 

0 

Variable 

in NUM 

Option 

exerci sed 


MTOTMX 

= 

0 

Striated flow solution. 


MTOTMX 

= 

1 

Fully mixed flow solution 

i (only for NCEC71 = 0 or 1) 

MODECH 

= 

1 

Input 

zone mass flows . 


MODECH 

= 

2 

Input ; 

zone fuel and oxidizer reservoir pressures. 

MODECH 

= 

3 

Input chamber static pressure. 

NAORW 

= 

1 

Input injector areas. 


NAORW 

= 

2 

Input injector velocities 

, (only for M0DECH=1 or 3). 

MGORLI 

= 

1 

Input injector gas constants. 

MGORLI 

= 

2 

Input injector densities. 


NSPPG 

= 

0 

Complete chamber solution 

1. 

NSPPG 

= 

1 

Special 

chamber solution 

(only for NCHA1=1, NCHA2-0) 

NCEC71 

= 

0 

Input thermochemical data by namelist EQU. 

NCEC71 

— 

1 

Calculates thermochemical 

1 data for frozen, calorical 




perfect flow. 


NCEC71 

= 

2 

Calculates thermochemical 

1 data for equilibrium flow. 

NCEC71 

— 

3 

Calculates thermochemical 

1 data for frozen flow. 
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Namelist EQU Variables (used only when NCEC71 = 0); 


Variable 

Dimension 

Type 

Common label 

Value before read 

T0G a 

20 

R 

THR 

by first-initial case 
undefined 

GAMG 

20 

R 

TWO 

undefined 

MWG 

20 

R 

THR 

undefined 

NCHANG 

1 

I 

— 

0 


NCHANG must be read for each case under namelist EQU when NCEC71=0. If 
an initial case is being considered NCHANG must be greater than zero and 
thermochemical data variables for each flow zone starting from the inner- 
most flow zone must be input under namelist EQU. For stacked cases when 
NCEC71=0 and new thermochemical data are not to change from the preceding 
initial case, only NCHANG=0 need be read under namelist EQU. For stacked 
cases when NCEC71=0, control card LETOUT may be input with any integer 
greater than one in columns 1 to 5 and parameters in namelist NUM will not 
be output by a namelist write. ^ 


a TzeroG 

^Namelist NUM print out for cases after the first case will not be completely 
applicable to the cases since variable parameters in NUM from the previous 
case will be seen. 
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Thermo Data: 




Card Order 

Contents 

Format 

Card Column 

1 

THERMO 

3A4 

1 to 6 

2 

Temperature range 

3F10.3 

1 to 30 


for 2 sets of coef- 
ficients: lowest, 

common, and highest 
T, # K 


3 


4 


5 


6 


Species name, Date 
Atomic symbols and 
formula 

Phase of species 
(S,L, or G for solid, 
liquid, or gas respec 
tively) 

Temperature range 
Integer 1 


3A4, 2 A3 
4(A2,F3.0) 

A1 


2F10.3 

115 


1 to 12, 19 to 24 
25 to 44 

45 


46 to 65 
80 


Coefficients a. 5(E15.8) 1 to 75 

(1=1 to 5) 1 

(for upper temperature 

interval ) 

Integer 2 15 80 

Coefficients a. 5(E15.8) 1 to 75 

(1=6 and 7) 1 

(for upper temperature 

interval) and a. (i=l to 3) 

(for lower temperature 
interval ) 

Integer 3 15 80 

Coefficients a. 4( El 5. 8) 1 to 60 

(1=4 to 7) 1 

(for lower temperature 

interval ) 

Integer 4 120 80 


(a) repeat cards 3 to 6 for each different species 

(Final END 3A4 1 to 3 

Card) 
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Reactants Data: 


Order Contents 


Format 


Card Columns 


1 REACTANTS 


3A4 


1 to 9 


2 


3 

4 


One card for each reac- 
tant species (maximum 15). 
Each card contains: 

(1) Atomic symbols and 

formula numbers (maximum 
of 5 sets) 9 . 

(2) Relative weight 

(3) Enthalpy , cal/mole 

(4) State: ,L, or G for 

solid, liquid, or gas, 
respectively 

(5) Temperature associated 
with enthalpy in (3),°K 

(6) F if fuel or 0 if oxi- 
dant 2 

(7) Density in g/crn (op- 
tional ) 

blank 

NAMELISTS 


5(A2 ,F7.5) 

1 to 45 

F7.5 

46 to 52 

F9.5 

54 to 62 

A1 

63 

F8.5 

64 to 71 

A1 

72 

F8.5 

73 to 80 

3A4 

1 to 9 


(c) repeat cards 1 to 4 

for each flow zone 


a Program will calculate the enthalpy (3) for species in THERMO data at 
tne temperature (5) if zeros are punched in card columns 37 and 38. 

b Relative weight of fuel in total fuel or oxidant in total oxidant. 
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APPENDIX F 

SAMPLE CASES AND SAMPLE INPUTS 


SAMPLE INPUTS: 


RESERVOIR PRESS. FIXED — I ZONE — FROZEN — N2H4/N204— PRF=30KPSFA ,PRO=30KPSFA 
LIQUID INJECTION — RC/RT=1.0 — CORE 0 /F= 1 . 5~NCEC7I = 0 

1 

CNUM PRMT=0.,1., .002, .001,0. ,1., RT=.l, RCRT=1., EPSC=6., ESPG=.50, 
OXFULE=l.5»19*0., PRF=30000. , 19*0. , PR0=30000. , 19*0 . , AITAC=.0005, 

AIPCT = l. ,19*0., RHOIFU=i. 9488, 19*0., T IFU= 537. , 1 9*0. , 

RHOI OX ~2. 7787 ,19*0. ,TI0X=537. ,19*0. , CD I NJ= . 7 , 19*0. , 

HTOTMX=0 *M0DECH=2 ,N AORW = 1 » MGORL I=2*NCHA1=1,NCHA2=0,NCEC71=0* CENO 
CEQU NCHANG= 1 , T 0G=5400 • ,19*0. , MWG=2l.l, 19*0. , G AMG= 1.23,19*0., LEND 

ZONE MASS FLOWS FIXED — 2 ZONE — FROZ.FROZ — MOTI (1) =. I , MDT I ( 2 ) = . 02 5 
GAS INJECTION — RC/RT=.5 — STRIATED FLOW — INJ VELOCITIES FIXED — 80/20 STRIATION 
2 

CNUM PRMT=0.,1.,.004, .001,0. ,1., RT=.l, RCRT=.5, MDT I =. 1 , .02 5, 1 8*0 . , 
0XFULE=6. ,.5,18*0., R IFU=2* 24665. , 18*0. , T I FU= 2*537 ., 18*0. , 

WIFU=170. ,600. ,18*0., R 1 0X = 2* 1 554. , 18*0. , T IOX= 2*53 7. , 1 8*0. , EPSC=4., 

WI 0X=1 70. ,600. ,18*0. , AG AMGU=1. 2092,1.3792, 18*0. , AT0GU=6089. , 1460. ,18*0., 
AMWGU= 12. 9 69, 3. 024, 18*0. , M0DECH=1 , MGORL 1=1 ,NCHA 1=1 ,NCHA2= 1, 

MT OT MX=0 « N A0RW=2 , NCEC71 = 0, LEND 
CEQU NCHANG=1 , T0G=6089 ., 1460 ., 18*0. , MWG= 1 2. 969 , 3. 024, 1 8*0. , 

GAMG=1. 2092, 1.3792, 18*0., CEND 

CHAMBER STATIC PRESS. FIXED — 2 ZONE — EQULIBRIUM — H2/02 AT 50 PSIA 
GAS INJECTION — RC/RT=.5 — CORE 0/F=6., SHEATH 0/F=.5— 80/20 STRIATION 
2 

CNUM PRMT =0. ,1.,. 008,. 00 1,0. , 1 . , RT= . 2, RCR T= . 5 , EP SC=4 . , PG = 7200. , 

T I F U=2 *53 7. ,18*0. ,R IFU= 2* 24665. , 1 0*0., TIOX= 2* 53 7. , 1 8*0. , 

R! 0X= 2*1 554. , 18*0. , W I FU= 1 50 . , 600. , 1 8*0. , W I 0X= 1 50. , 600. , 18*0. , 

PCTM0=.8,.2, 18*0. ,MT0TMX=0, MODE CH=3,N AORW= 2, MGORL 1=1, NCHA 1=1, NCHA 2=1, 


NCEC7l=2fOXFULE=6* , .5, 18*0. » 

L END 
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GAS INJECTION — RC/RT=. 5— CORE 0/F=6., SHEATH 0 /F= . 5--60/40 STRIATION 
71 

CNUM PCTM0=.6 ,.4, 18*0. , CEND 
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CHAMBER STATIC PRESS. FIXED— 2 ZONE— EQUL I BRIUM— H2/02 AT 50 PSIA 
CAS INJECTION— RC/RT-. 5 — CORE O/F-6.# SHAATH 0/F-.5— 90/20 STRIAT10N 
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PRSPFD 


APPENDIX G 
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Subprograms added to CEC71 


APPENDIX G (Continued) 
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ROUTINE TREE DIAGRAM - MODULE II | SAVE 
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Flow diagram for ma i n routine 
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Flow diagram for subroutine CHAMBR 
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TABLE I - PROGRAM INPUT FOR 


PRESELECTED THERMOCHEMICAL 
INPUT DATA 


Comment 

Card 


Comment 

Card 


Integer 

Control 


Card (LETOUT) 


✓ — — - 

71 

' & NUM 

Geometry data 
Integration data 
Chamber operation mode data 
or mass generation data 
Option data 
Injector data 
Initial estimations data 

for variable chamber pressure 

V 



190 






TABLE II - PROGRAM INPUT FOR 
INTERNALLY GENERATED THERMOCHEMICAL DATA 


— 

Comment 

Card 
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Comment 

Card 


Interger 

Control 


Card (LETOUT) 


& NUM 

Geometry data 
Intergration data 
Chamber operation mode data 
Option data 

Injector data 


THERMO 
Code Card 


THERMO 
data cards 


OR 



THERMO 
data on 
disk 4 


OR 


THERMO 
data on] 
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Code Card 
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ONE SET REQUIRED FOR EACH ZONE. 
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Figure la. Schematic of Nozzle Chamber. 
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Figure lb. Schematic of Nozzle Throat. 
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Figure 3. Discharge Coefficient Comparison Between 
Experimental and Analytical Results. 
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Figure 4. Throat Plane Mach Number Distribution 
Comparison Between Experimental and 
Analytical Results. 
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Figure 5. Throat Plane Pressure Ratio Distribution 
Comparison Between Experimental and Anal- 
ytical Results. 
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Figure 6. Variation of Discharge Coefficient with 
Cooling Mass Ratio for R c /R-p = 0.5. 
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Figure 7. Variation of Discharge Coefficient with 
Cooling Mass Ratio for R c /Rj = 2.0. 
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Figure 10. Discharge Coefficient for Low R c /Ry with 
Various Amounts of Cooling Mass Flow. 
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Figure 12. Characteristic Velocity For Mixed and Unmixed 
Flow With Chamber Static Pressure Fixed. 
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Figure 13. Characteristic Velocity For Mixed and Unmixed 
Flow With Reservoir Pressures Fixed. 
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